2D Backward Stochastic Navier-Stokes Equations 
with Nonlinear Forcing * 

Jinniao Qiu \ Shanjian Tang ^\ and Yuncheng You ^ 

January 20, 2013 



Abstract 

The paper is concerned with the existence and uniqueness of a strong so- 
lution to a two-dimensional backward stochastic Navier-Stokes equation with 
nonlinear forcing, driven by a Brownian motion. We use the spectral ap- 
proximation and the truncation and variational techniques. The methodology 
features an interactive analysis on basis of the regularity of the deterministic 
Navier-Stokes dynamics and the stochastic properties of the Ito-type diffusion 
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1 Introduction 



Let {fl, ^ , {^t}, P) be a complete, filtrated probability space, on which defined 
is a standard 1- dimensional Brownian motion {Vri}j>o, whose natural augmented 
filtration is denoted by {^i,t G [0,T]}. We denote by CP the cr-Algebra of the 
predictable sets on x [0,T] associated with {^t}t>o- The expectation will be 
exclusively denoted by E and the conditional expectation on will be denoted by 
E^^^. We use a.s. to denote that an equality or inequality holds almost surely with 
respect to the probability measure P. 

The theory of backward stochastic differential equations (BSDEs) has received 
an extensive studies in the last two decades in connection with a wide range of 
applications as in stochastic control theory, econometrics, mathematical finance, and 
nonlinear partial differential equations. See p]-[H[71[8| [m[22] for details. In Tang [T8] , 
a very general system of backward stochastic partial differential equations (BSPDEs) 
are studied. However, they are semi-linear, and thus exclude the nonlinearity of the 
well-known Navier-Stokes operator. In this paper, we concentrate our attentions to 
study the backward stochastic Navier-Stokes equation (BSNSE). 

The standard deterministic Navier-Stokes equation describing the velocity field 
of an incompressible, viscous fiuid motion in a domain of M."^ [d = 2 or 3) takes the 
following form: 

r dtu - uAu + (u ■ V)u + Vp + / = 0, t > 0; , . 

\ V ■u = 0, u{0) = Mo, ^ ' ^ 

where u = u{t,x) represents the (i- dimensional velocity field of a fiuid, p = p(t,x) is 
the pressure, z/ G (0, oo) is the viscosity coefficient, and / = f(t,x) stands for the 
external force. Let {u,p) solve the equation (11. ip . By reversing the time and defining 

u(t,x) = ~u(T — t,x), p(t,x) = p{T — t,x), for t < T, 

then {u,p) satisfies the following backward Navier-Stokes equation: 

j dtU + uAu + {u-V)u + Vp+ f = 0, t <T; 

\V-u = 0, u{T)=uo. ^^-^^ 

Note that the time-reversing makes the original initial value problem of (11. ip become 
a terminal value problem of (II. 2p . 

We shall study the following two-dimensional backward stochastic Navier-Stokes 
equations (briefiy 2D BSNSE) in with a spatially periodic condition and a given 
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terminal condition at time T > 0: 

' du{t, x) + {uAu{t, x) + {u ■ V)u{t, x) + {a- V)Z{t, x) + Vp{t, x)} dt 

= -f{t, X, u, Z) dt+Z{t, x) dWt, {t, x) G [0, T) x R^; 
< divu{t,x) = 0, G [0,T) X R2; (1.3) 

u{t,x + aci) = u{t,x), G [0,T) X i = l,2; 

u{T,x) = ^{x), {x,u) G M^. 

Here a is a measure of "correlation" between the Laplace and the Brownian motion, 
{61,62} is the canonical basis of R^, a > is the period in the ith direction, u = 
(mi,M2) is the random two-dimensional velocity field of a fluid in M^, / represents 
the external forces which allow for feedback involving the velocity field u and the 
stochastic process Z and may be inhomogeneous in time. The terminal status of 
the velocity field is a known random field ^ on the underlying probability space. 
For notational convenience, however, the variable a; G in various functions and 
solutions will often be omitted. 

It is worth noting that, though sharing the same name, our BSNSE essentially 
differs from that of Sundar and Yin p^7] since the sign of the viscous term "z/Am" 
differs. Furthermore, we allow the external force / to depend on both unknown 
fields u and Z in a nonlinear way, and the drift term to depend on the gradient of 
the second unknown field Z. 

In |9l|10], Cauchy problems for the (forward) stochastic Navier-Stokes equations 
in R*^ driven by a random nonlinear force and a white noise are studied and the 
existence and uniqueness of a global martingale solution have been proved. As a mo- 
tivation BSNSEs emerge in regard to inverse problems to determine the stochastic 
noise coefficients from the terminal velocity field as observed. In [6l[23], a stochastic 
representation in terms of Lagrangian paths for the backward incompressible Navier- 
Stokes equations without forcing is shown and used to prove the local existence of 
solutions in weighted Sobolev spaces and the global existence results in two dimen- 
sions or with a large viscosity. In [ISIITT], the existence and uniqueness of adapted 
solutions are given to the backward stochastic Lorenz equations and to the backward 
stochastic Navier-Stokes equations (11.31) in a bounded domain with a = 0, z/ < and 
the external force f(t,y,z) = f(t). 

The rest of the paper is organized as follows. In Section 2, we introduce some 
notations, assumptions, and preliminary lemmas, and state the main result (see 
Theorem 2.1). In Section 3, we consider the spectral approximations and give relevant 
estimates. In Section 4, we prove the existence of an adapted solution to the projected 
finite dimensional systems for our 2D BSNSE. Finally, in Section 5, we give the proof 
of Theorem 2.1. 
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2 Preliminaries and the main results 

Let G = (0, a)^ be the rectangular of the period. For any nonnegative integer m, we 
denote by H"^{G) the Sobolev space of functions which are in L^{G), together with 
all their derivatives of orders up to m and by H^{G) the space of functions which 
belong to HJ^^iM?) (i.e., u\o G if™(0) for every open bounded set 0) and which are 
periodic with period G. H^{G) is a Hilbert space for the scalar product and the 
norm 

{U,v)m= ^ / D°'u{x)D°'v{x) dx, \u\m = {{U,u)my^'^, 

where a = (ai, 0:2), [a] = ai + 0:2, and 



dxi^ dx 
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with ai and ^2 being two nonnegative integers. The elements of HJ^{G) are charac- 
terized by their Fourier series expansion: 

H;', = {u: u=J2 Cfee2-'=-/^ = c_fc, ^ \k\'^'\ck\' < 00}, (2.1) 

1 /2 

and the norm is equivalent to the norm 

{Efcez2(l + \k\''"')\ck\T . Set a Hilbert 

subspace of H^{G): 

H^^{G) = {ue H^e{G) : in its Fourier expansion ([2l]), cq = 0}, (2.2) 

with the norm \u\m,o = {J^kez^ \k\'^"'\ck\'^V^'^- Actually, through ( 12. ip and ( 12. 2p . we 
can define Hp^{G) and H^^iG) for arbitrary m G M. Moreover, H^^{G) and H-^iG) 
are in duality for all m G M. 

As in the framework of treating the deterministic Navier-Stokes equations (c.f. 
[UlEniEI]), we set up three phase spaces of functions of the spatial variable x & G 
as follows: 

H = {ife Hl^{G) X Hl^{G) : divyp = in M'}, 
V = {^e Hl^{G) X Hl^{G) : div(p = in M^}^ 

r = i/p|(:7~(R2) X c°°(r2). 
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where C°°(]R^) is the set of smooth functions on M^. Then both H and V are Hilbert 
spaces equipped with the respective scalar product and norm 

2 

(0, : = UWh := {(0, cI>)hY'\ <P,^eH; 

(2.3) 

(0, : = X^pj0\ /^,V)o, UWv ■■= {(0, 0,¥^ e 

For simphcity, we denote || ■ || and (■, ■)h by || ■ || and (■, ■) respectively. The dual 
product of -0 G and ip G V will be denoted by (-0, f)v',v and it follows that 

(</>i, (f)2)v',v = (f)2) 01 e V, 02 e if. 

For a little notational abuse we still denote by (■, ■) the dual product (■, ■)v',v- We 
shall use | ■ | to denote the absolute value or the Euclidean norm of . The set of all 
positive integers will be denoted by Z"*" or N. The Lebesgue measure of the domain 
G will be denoted by IG]. Define := H^^KG) x if^(G). 

For any finite dimensional vector space F and a,b & F, we denote by a ■ 6 the 
scalar product on F. Throughout this paper, we assume that the external force term 
f{t, u, Z) and the terminal value term ^ are EI°-valued and H^-valued respectively, so 
the solution pair [u, Z) of (11.31) must be EP x H°-valued. By applying the projection 
^ H (see [21]), since 

= {ifj eif -.tp = Vq for some q E HI^{G)}, 

we can formulate the above terminal value problem of the 2D BSNSE (ll.3p into 
the following problem to solve the backward stochastic evolutionary Navier-Stokes 
equation, 

-du{t) ={-uAu{t) + B{u{t)) + JZ{t) + /(t, u(t), Z{t))}dt 

-Z{t)dWt, tG[0,T), (2.4) 

where 

U{u, v) = F{{u ■ V)v) -.VxV ^V, B{u) = U{u, u):V V, 

2 

JZ = F{{a-V)Z), a{t,x) = {a\t,x), a\t,x)), ((a ■ V)Z)^= a^Z^^, 
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and 

whose domain is D{ A) = M.^ f]H and by Poincare inequality we can show that V — 
D{A^/'^). Accordingly we shall adopt the equivalent norm \\(p\\v — ||V</?|| = ||^^/^(/?||. 
Then all H, V, and D{A) (with the graph norm || • ||d(A)) are separable Hilbert spaces. 
With a little notational abuse we still use / and Z for the projections P(/) and P(2'), 
respectively, 

Note that the above Stokes operator A : D{A) H is positive definite, self- 
adjoint, and hnear, and its resolvent is compact. Therefore, all the eigenvalues of 
A can be ordered into the increasing sequence {Ai}^^. The corresponding eigen- 
functions {ej}^^ form a complete orthonormal basis for the space H, which is also 
a complete orthogonal (but not orthonormal) basis of the space V. With the iden- 
tification H — H' by the Riesz mapping, one has the triplet structure of compact 
(consequently continuous) embedding, 

V CH cV. (2.5) 

In what follows, C > is a constant which may vary from line to line and we 
denote by C(ai, a2,- ■ ■) or Cai,a2,--- a constant to depend on the parameters Oi, 02, • • • . 
We make the following three assumptions. 

Assumption (Al). The i?- valued mapping / is defined on Q x [0,T] x \^ x 
H and for any {u,z) e V x H, f{-,u,z) is a predictable and if-valued process. 
Moreover, there exist a nonnegative constant ^ and a nonnegative adapted process 
g G L°°((], L^([0, T])) such that the following conditions hold for all v,vi,V2 G V, 
0,01,02 G H and (w,t) G x [0,T]: 

(1) . the map s h- )■ {fit, vi + 5^2, 0), v) is continuous on R; 

(2) . 

{fit, Vi, 0i) - fit, V2, 02), Vi - V2) 
< piv2) {\\Vi - V2IP + \\Vi - ^211(1101 - 02|| + ll^^l - ^^2||v)) ; 

where p : y — >■ (0, -|-oo) is measurable and locally bounded; 

(3) . 

{fit, V, (t>),v)< git) + 6||0||2 + Qie)\\vf + 

where ^ : (0, 1] ^ M+ is continuous and decreasing; 

(4) . 

\\fit,v,ci>)f<{git)+mv\\i+ur))piiv). 

where pi : H ^ (0, -|-oo) is measurable and locally bounded. 
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Remark 1. In fact, (1) and (2) of Assumption (Al) implies f{t,x,u,z) is locally 
Lipschitz continuous with respect to z in the following sense: 

\\f{t,u,z) - f{t,u,Z)\\v^ < C{\\u\\v)\\z - Z\l 

for all {u, t) enx[0, T],u eV and z,Z e H. Actually, for any (peV- {0}, e e R+, 

{fit, u + 60, z) - fit, u, Z),ct>)< Ci\\u\\v)mv{em + \\z - Z\\}. 

Letting e J, 0, from the arbitrariness of we conclude that the local Lipschitz conti- 
nuity holds. 

Assumption (A2). The function defined on x [0, T] is real- valued CP- 
measurable such that Id-' I < A, almost surely for j' = 1,2 and all t e [0,T], for 
some A e (0, oo). 

Assumption (A3) (super-parabolicity). There exist two constants A > and 
A > 1 such that 

2Her-P(<7(t)-0'>2A|er 

holds almost surely for all (i,^) e [0,T] x M?. 

Note that, in Assumptions (A2) and (A3), our cr is defined independent of the 

spatial variable x. 

For Banach space B and p > 1, define 

jSf|(0, T; B) := {0 e i7(Q x [0, T];B)\ {0(-, i)}o<t<ris a predictable process}. 
Define 

M[0, T] := L%iQ; C([0, T]; if)) n JSf|(0, T; F) 

and 

:= ^[0, T] := M[0, T] x jSfJ:(0, T; H) 

equipped with the norm 



















sup ||M(t)|P 




/ \Ht)fvd't 


+ E 


Hmwdt 


) 




_te[o,r] 




Jo 




Jo 





1/2 



Throughout the paper, define 

^t, (f>, cp) := -uA<f> + Bi<f>) + Jip + fit, (f>, cp), (0, ip)eVxH. (2.6) 
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Definition 1. (weak solutions) For C, € L^^{Q] H) given, we say tliat {u,Z) G ^ 
is a weak solution to (12 .4^ if for any ip & there holds almost surely 

(M(t), ip) = V)+ [ ms,u{s),Z{s)), if) ds 

(2.7) 



^ {Z{s), if) dWs, Vte [0,T]. 



Definition 2. (strong solutions) For ^ G L'i^^{Q;V) given, we say that {u,Z) is a 
strong solution to (12.41) if {u, Z) is a weak solution and 

(u, Z) G (L2^(^]; C([0, T]; r)) n ^|-(0, T; D{A))) x ^|(0, T; V^). 

Remark 2. If we have verified that (m, Z) G (L^(fi; C([0, T]; 1/)) n^|(0, T; D(A))) x 
^1(0, T;1/) and that 

u{t)=i + <!){s,u{s),Z{s)) ds- Z{s) dWs a.s. in H, 

By the stochastic Fubini theorem (see, [T2l Theorem 4.18]), we can check that (u, Z) 
is a strong solution to (12.41) . 

The main result of the paper is stated in the following theorem. 

Theorem 2.1. Under Assumptions (Al)-(A3), for any ^ G Ljp^(fi;V), the 2D 
BSNSE problem (12. 4p admits a unique strong solution such that 



esssup [[^(s)!!^ + 

{iv,s)enx[o,T] 



(2. 



< |ll5'l|L-'(n,Li([o,T])) + UWl^^in-y) 

where C is a constant depending on z/. A, A, /3, g, pi and T . 
For the trilinear mapping 

2 2 „ 

b(u, V, w) := (n('U, v), w) = y > Ui -—-^ Wj dx, u,v,w G H. 

i=i j=i * 

we have the following instrumental regularity properties. 
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Lemma 2.2. The following properties hold for any two-dimensional bounded domain 
G, where Cq is used to denote different constants only depending on G. 

|,/ M ^ ol/2|| ||1/2|| ||l/2 1/2,, ,,1/2 ^ -J. 

\b{u,v,w)\ < CG\\u\\]i^\\Au\\]i^\\v\\v\\w\\H, u e D{A),v e V,w G H, 

\b{u,v,w)\ < GG\\u\\f\\u\\\!'\\v\\\!'\\Av\\f\\w\\H, ueV,ve D{A),we H, 

1 /2 1/2 

\b{u,v,w)\ < C'G||M||_ff||^||y II^^IIh ! u £ H,v G G -D(y4). 

Moreover, 



(n(u,f), w) = — (n(u, ty), t>), forM,i;,wGV, 
{Il{u,v), v)=0, foTU,VGV. 



(2.10) 



For u G -D(^), tfe /iawe B{u) G if, 

< CG\\u\\f\\u\\v\\Au\\]l\ (2.11) 

anc? especially, for the periodic case, we have 

{Il{v,v), Av) = 0, for V G D{A) (c.f JM Lemma 3.1, Page 19]). 

The proof of Lemma 12.21 can be found in [T9l|20]. The first inequahty in (12. 9p 
will be most useful in this work and its coefficient equals 2^/^ for any bounded and 
locally smooth domain in space dimension n = 2, which was proved in [191 Lemma 
3.4]. The following lemma shows the regularity of functions in Hq{G) for a 2D 
domain G, whose proof is available in [19]. 

Lemma 2.3. For any two-dimensional open set G, we have 

\\v\\lHg) < ^'^'\\v\\%) W'^^WlHgv ^ ^ ^oiG). (2.12) 

We have the following two versions of Gronwall- Bellman inequalites, whose proof 
is referred to [2l[5]. 

(The Gronwall-Bellman Inequality): If a nonnegative scalar function g{t) is con- 
tinuous on [0, T] and satisfies 

9{t)<{>)g{T)+ j {ag{s) + h{s))ds, tG[0,T], (2.13) 
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where a > is a constant and h : [0, T] — )■ M is integrable, then 

9{t) < (>) e"(^-*)^(T) + e^'^'-'^his) ds, t G [0, T]. (2.14) 

(The Stochastic Gronwall-Bellman Inequahty): Let {fl, 5", F, P) be a fihered prob- 
abihty space whose filtration F = {3^t '■ t G [0,T]} satisfies the usual conditions. 
Suppose {Yg} and {Xg} are optional integrable processes and a is a nonnegative 
constant. If for all t, the map s ^ E[Ys\3^t] is continuous almost surely and 



Yt < {>)E 
then we have almost surely 

Yt < (>)e"(^-*)E[yT|3't] + E 



(X, + aY,) ds + Yt 



e"('-*)X, ds 



Vt e [o,T]. 



In this paper, we prove the existence and uniqueness of an adapted solution to the 
terminal value problem f ll.3p of a two-dimensional backward stochastic Navier-Stokes 
equation with nonlinear forcing and the random perturbation driven by the Brownian 
motion. We use the spectral approximation, combined with the truncation and 
variational techniques, which is also a kind of compactness method. The methodology 
features an interactive analysis based on the regularity of the deterministic Navier- 
Stokes dynamics and the stochastic properties of the Ito-type diffusion processes. 



3 Spectral approximations and estimates 

In this section we consider the spectral approximation of the problem (12. 4 p obtained 
by orthogonally projecting the equation and the terminal data on the finite dimen- 
sional space 

Hn = Span {ei, 62, ■ • • , cat}. 

Define 

N 

Pn-.V ^ Hn, PnI = J2{f^ e,)e„ / e V . 

1=1 

Then ||Pjv/|p = is the orthogonal projection on Hjy, which is 

called the spectral projection. It is worth noting that || ■ || and || ■ ||y are equivalent 
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in Hn and that Hn = Vn '■= PnV. Define 

N 

= PnA, B'^iu) = PnB{u), J^'Z = PmJZ := ^(JZ,e,)e,; 

r (■ ■ ■ ) = PNf{- ■■), Z'^it) = P^Z{t), and e = PnC 

Then the projected, A^- dimensional problem of approximation to the problem (12 .4^ 
is defined to be 

+ Z''{t)dWt, tG [0,T); 

(3.2) 

Note that the projection does not affect the Brownian motion {Vrt}t>o, and also 
that, the finite dimensional approximation equation (13. 2p does not satisfy the condi- 
tions listed in p!]. 

We shall conduct a priori estimates for the adapted solution to the finite dimen- 
sional approximation problem (13.21) . 
First, by means of Young's inequality 

ab < -aPeP + — 6^ 1 + 1 = 1, ab>0,e>0, 
p qe'i p q 

under the Assumptions (Al)-(A3), we have 

2($(t,0,yp), 0)-||^f 
= 2{-vA<P + B{<P) + Jv + fit, 4>, <^), </)) - llv^f 
= - 2u{A<P, 0) - 2(/(t, 0, y^), 0) - 2{^, (a ■ V)0) - yf 
< - 2uUry + 2ig{t) + eM' + Q^H' + PUIWI 



+ 2^\\ip\\ \\\{cr ■ V)0|| — ||<y3||^ (choose e small enoug h) (3-3) 
A 



— \\ipf + 2q(t) H ^ 

|2 -'-||,,l|2 I n„/4.\ , r^llJ.ll2 



< - 2Aii0f^ - ~^yr + m) + + aii^ii^ + cm 



= - Mmv - -jy^w^w + m) + cm\ (0, v^) g v x f^, 

where the constant C depends only on A, A, g and /3. 
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Lemma 3.1. Let the conditions of Theorem 12.11 hold. If {u {■) , Z {■)) G 
adapted solution of the problem i\3.2L then we have almost surely 
(!)■ 



IS an 



sup \ \\u 

te[o,T] 



N 



u'^isWv+WZ^'isW ds 



N I 



(3.4) 



< C (||fi'||L°°(C,Li([0,T])) + K\\L<^{n,H)) ) 

where C is a constant depending only on T, u, A, A, /3 and g; 
(2). 



sup < \\u^{t)\\l + E^ 



te[o,T] 



\\Au''{s)r+\\z\s)rvds 



(3.5) 



\V ^ ^^t 
= (n,Li([0,T])) + ll^llLoo(Qy)) 

with C being a constant depending only on u, A, A, /3, g, pi and T . 

Proof. Applying the backward Ito formula to the scalar-valued, stochastic process 
||M^(t)|p, and noting that 



(5^(n^(t)),n^^(t))=0 



N I 



we have 



lie 



2u {A'^u'\s),u'\s))ds + 2 {f'\s,u'\s),Z'\s)),u'\s))ds 



+ 2 l{j''Z^{s),u'\s))ds-2 l{Z'\s),u'\s))dW, 



.N 



.N I 



Z^isWds 



lie 



2u {Au'\s),u'\s))ds + 2 {f{s,u'\s),Z'\s)),u'\s))ds 



+ 2 l{JZ''{s),u'\s))ds-2 l{Z'\s),u'\s))dW, 



.N 



.N I 



Z^'is^fds. 
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In view of f l3.3p . we have 



N\\2 



(2($(.,«^(.),Z^(.)), «^(.))-||Z^(.)f) ds 



-2 / {Z'\s), u''{s))dW, 



+ 



T 



X2 _ T 



Z^^(.)f + 2^7(.) + C||n^(s)f ) ds 



4A2 



where the constant C is independent of A^. Since 



E 



sup 



< 2E 



< CE 



< CE 



sup 

re[i,T] 



u'^is)) dW^ 
{Z^{s), u^{s)) dW, 



T \ 1/2- 

WZ"^ is)\W is)\? ds 



sup 

sg[i,T] 



T \ 1/2- 











< (l/2)i? 


sup ||«^(s)f 


+ CE 


j^\\z^{s)rds 




sG[t,T] 







(3.6) 



(by BDG inequahty) (3.7) 



taking the conditional expectation on both sides of the second inequahty of (13.61) . 
we obtain 



\u 



N(+\\\2 



\E. 

<E,^,\\if\^CE,^, 



\n''{s)fyds 



4A2 

W)^\\u''{s)f) ds 



Z^(s)||2 ds 



a.s.. 
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Prom the stochastic Gronwall- Bellman inequality, it follows that 



^""{8)111(18 



\\z^{s)rds 



sup \\\u''{t)\\' + Ej, 

< C (||5'||L-(n,Li([o,T])) + 

< C (||5'||L-(n,Li([o,r])) + ||Clli°c(n,H)) , a.s. 



where C is a constant depending only on T, i/, A, A, /3 and g. 

On the other hand, as (-B(u), Au) — 0, using Ito formula, we have 



(3.8) 



\u 



N 



Wily = 11^^11^-2^ 



T 



+ 2 / (J^Z^(s), A"u''{s))ds 



N„.Ni 



.N 



sn'yds 



.N, 



-i: 



WZ'^isWvds 



rT 2 

-'I ^ 

i=i 



(V(Z^r(.), {a-V)V{u''ns)))ds 



2 j\z''{s), Au''{s))dWs, t e [0,T]. 



By the integration-by-parts formula and the fact that the integrals on the boundary 
dG of G vanish, we obtain 



ll^0f = E llV'^ilv> V0eD(A). 



i=l 
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It follows that 



-2v I \\Au''{s)rds-2 I ^(V(Z^r(5), {a-V)V{u'')\s)))ds 
Jt Jt 

2 T 2 

< -2uJ2 [ \\Viu''is)y\\lds+ [ ^||A(a-V)V(«^r(.)|prf. 

i=i -^t -^t i=i 

< -2\Y,j\v{u''{s)nlds + \-^j\z''{s)\\lds, a.5.,Vte[0,r]. 



Therefore, we have 



-2A \\Au''{s)fds + 2 j'^{f{s,u^{s),Z^{s)), Au"" {s)) ds 



<llell^-2j[ (Z'^is), Au''is))dW,-^J^ WZ^'isWyds 

- X WAu'^is) fds 11/(5, u^'isl Z^s))f ds 

< Urv-2l\z^{s), Au^{s))dWs-^^ \\Z''{s)\\lds 



1 '■^ 



-A \\Au'^{s)rds + -J^ [g{s) + mu''{s)rv + \\Z''{s)r)]p^{u'')ds 
< m\l-2j\z'^{sl Au'^{s))dWs-^l^^\\Z'^{s)\\lds 

-A r WAu'^isn'ds + C r [gis)+m^^''is)\\l + \\Z''is)r)] ds. 
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As 



E 



sup 

re[t,T] 



< 2E 



< CE 



sup 

re[t,T] 



{Z''{s),Au^{s))dW, 

E [\^{z''y{s),v{u'^y{s))dWs 

1/2- 



1=1 



(by BDG inequality) (3.10) 



sup \\u^{s)fy 

se[t,T] 



CE 



< C{N) |e 


sup ||u^(s)||y 


+ E 




s&[t,T] 





taking conditional expectation on both sides of (13.91) . we get 



\\u''{t)fy + XE.^, 

<E,^^[Ufv]+CE^, 



\\Au'' {s)\\' ds 

T 



+ 



2A2 



-E 



{9{S) + Wu"^ {sWy + \\Z^\s)r) ds 



N I 



In view of (13.81) . we conclude that, with probability 1, 



sup <^ \\u 



N 



V + Ej^ 



\\Au^{s)fds 



IZ'^isWyds 



(3.11) 



(3.12) 



< C (||5'||L°°(f7,Li([0,T])) + ll^llL-'(Q,y)) , 

where C is a constant depending only on z/, A, A, /3, g, pi and T. □ 
Lemma 3.2. For any u,v & V and (j),ip E H , 

\{B{u)- B{v),u-v)\ < h\u-v\\l + h\v\\l\\u-v\\\ (3.13) 

4 A 

Moreover, under Assumptions (A1)-(A3), there exists a positive constant K depend- 
ing on A such that 



-2{^{t,u,(t))-^{t,v,^),w) 
+ ||wf (k+j\\v\\1 + Kp^{v] 



A^ + l 
2A2 ' 



— II 2 II I 

w\\ > —\\w\ 



(3.14) 
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holds almost surely for any t ^ [0,T], u,v & V and ip,(f) E H with w := u — v , w 
— ^9, and $ being defined by (12.61) . Define 

r,{t) = (^K + j\\u{s)fy + Kp''{u{s))^ ds 



and 



r2{t) = (K+j\\vis)\\l + Kp'ivis))^ ds, 



for arbitrary u,v E Lj^{Q; L'^{0,T;V)) , and let w{-) = u{-) — v{-). Then for any 
E L'^i^] -^^(0, T; H)) and w{-) := 0(-) — we have for i = 1,2 

- (2<l>(t, u, (f)) - 2$(t, V, if) + ^^w, w) + ii^ll^f > 0, a.s.. (3.15) 

Proof. Let w = u — v. Then 

{B{u) - B{v),u-v) 
= - {Il{u,w),u) + {n.{v,w),v) 
= - (n(M, w),v)^ (n(t;, w),v) = -{B{w), v). 

By the first inequahty in Lemma 12.21 we can get 

\{B{u) - B{v),u- v)v',v\ = \{B{w),v)\ = \{Il{w,v),w)v',v\ 

< 2^^^\\u — — ^||y||f||y < jll^i — ^||v' + -rll'"~'^ini'^llF- 

4 A 

It follows from the Assumptions (A1)-(A3) that 



- 2u{Aw, w) + 2{Jw, w) + 2{f{t, u, 0) - f{t, v, ^),w) 

< - 2u\\wfy + 2{w, (a ■ V)w) + 2p{v)\\wf + 2p{v)\\w\\{\\w\\v + \\w\\) 

3A.. 11 9 1.. 11 9 ,^ l-(-A^,|| 1,9 

< - yll^^llv + Jill^f + (1 - (3.16) 
+ {CC\)p'iv) + 2p{vmwr 

< - yikirv^ + ^ii^ir + + Kp^vmwf, 

where the constant K only depends on A. Hence, in view of fl3.13p . we obtain (I3.14p . 

Then (13.151) for i = 2 follows from (13.141) by direct calculation. The case of i = 1 
in (13.151) is shown in a similar way. □ 
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4 Solutions of the finite dimensional systems 

In this section, we consider the existence of an adapted solution to the projected, 
A^- dimensional problem fl3.2p of the 2D backward stochastic Navier-Stokes equations 
which we also call the finite dimensional system. To solve the finite dimensional 
system (13. 2p . we shall make use of the result of Briand et al. [1]. 

Consider the following backward stochastic differential equation (BSDE in short): 

Y{t) = C + g{s,Y{s),q{s))ds- q{s)dWs, (4.1) 

where C, is an R^-valued ^^^-measuable random vector and the random function 

g: [0, T] X ^] X X ^ 

is y X ^(R^) X =^(R^)-measurable. 

The following lemma comes from [T], Theorem 4.2]. 

Lemma 4.1. Assume that g and ( satisfy the following four conditions. 
( CI ). For some p > 1, we have 



E 



|Cr+ / |^7(t, 0,0)1 ds 







< oo. 



(C2). There exist constants a > and /i G R such that almost surely we have for 
each (t, y, y', z, z') G [0, T] x R^ x R^ x R^ x R^, 

\g{{t,y,z)-g{t,y,z')\<a\z-z'\, (4.2) 

{y-y\9ii,y,z)- g{t,y',z)) <iJ,\y-y'\^ (monotonicity condition) . (4.3) 

( C3). The function y ^ g(t, y, z) is continuous for any (t, z) G [0, T] x R^. 
(CJi^). For any r > 0, the random process 

Mt) ■■= sup \g{t,y,0)-g{t,0,0)\, t G [0,T] 

\y\<r 

lies in the space L^{Q x [0,T]). Then BSDE (14.11) admits a unique solution {Y,q) G 
LP(fi,C([0,T],R^)) X £P(^],L2([0,^],R^)). 

Remark 3. It is worth noting that our finite dimensional system does not satisfy the 
monotonicity condition (C2). In fact, by Lemma [3.21 our finite dimensional system 
only satisfies a local monotonicity condition in some sense, which prevents us to 
directly use this lemma to our finite dimensional system. 
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Lemma 4.2. For any M,N G Z^, define the function of truncation Rm{-) to be a 
function on Hn such that for X = XiCi, 

{1, if||X||<M; 
G(0,1), ifM< ||X|| <M + 1; 
0, if||X||>M + l. 

Thus Rm{') is uniformly Lipschitz continuous. For each n G denote (fn{z) = 
zn/{\\z\\ V n), z G and set 

71 

nM[t) V n 

where 

huit) = 4 \ {git) + I3Cn{M + l)^) esssup \pi{w)\ + C^,,(M + 1)^ 
> esssup ty, 0)11 

||'!«||<A/+l 

and Km G L^{Vt x [0,r]). Then under Assumptions (A1)-(A3), <|)^'*^'" satisfies the 
conditions (C2)-(C4) in Lemma W7\\ 

Proof. Under Assumptions (Al)-(A3) and Remark [1], we only need verify (14.31) . i.e., 
there is a uniform constant CN,M,n > such that 

($^'^'"(t,X,Z)-$^'*^'"(t,F,Z), X-Y)<C^,M,n\\X-Yf, a.s., (4.5) 

for any X,Y, Z & Hjy and all t G [0, T]. For any X,Y & H^, inequality (14. 5 p holds 
trivially if ||X|| > M + 1 and > M + 1. Thus, it is sufficient to consider the 




4) 
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case of < M + 1. We have 

X, Z) - F, Z), X-Y) 



Rm{X) 



n 



hM{t) V n 



($(t,X,(^„(Z))-<l>(t,F,^„(Z)), X-Y) 



n 



{Rm{X) - RM{Y)){^^''''''{t, Y, ifniZ)), X~Y) 



+ 



< [K+j\\Y\\^ + Kp\Y)]\\X-Y\ 



hM{t) V n 

( by (EUD of Lemma [32D 



(4.6) 



+ Cm\\X-Y\ 



4 



n 



hM{t) V n 



^''''''-{t,Y,^^{zm 



< (K+j-\\Y\\'^ + Kp^{Y)]\\X-Yf 



+ Cm\\X-Y\ 

< CAI^N,n\\X — Y\\^ 

which completes the proof. 



n 



hM{t) V n 



{huit) + Cn,m ■ n) 



□ 



Theorem 4.3. Let Assumptions (Al)-(A3) hold. For any C, € L'^{Q;V), the pro- 
jected problem (13. 2p admits a unique adapted solution {u^ {■), Z'^ {■)) G ^ for each 
given positive integer N, which satisfies 



<c{i+E[ur]} 

where C is a constant independent of N. 



(4.7) 



Proof. Step 1. Let us verify the uniqueness part. Suppose (m^, Z^) and (f ^, F^) 
are two solutions of the projected problem (13.21) . Note that the a priori estimates in 
Lemma [3.11 holds for both {u^,Z^) and {v^,Y'^). Denote by {U^,X^) the pair of 
processes {u^ 



- v^,Z^ -Y^). Define 



r{t) :-- 



ds. 
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An application of Ito formula and Lemma 13.21 yields 

eK*)||f/^^.Mi2 



\\X 
T 



4 



U^{s)r[K+-\\v-{s)ry + Kp\v-{s)) 



ds 



(4. 



< 



2A2 



Taking conditional expectations on both sides, we get 



+ 



2A2 



fi^)\\X''{s)\\^ds 



< 0, a.s., for any t G [0,T], 



which implies the uniqueness. 

Step 2. For any A^, M, n G Z+, following Lemma [4. 2 ^ we can verify that the pair 
(^^, $^'*^'") satisfies the conditions (C1)-(C4) in Lemma [4.11 Hence, by Lemma |4?T] 
there exists a unique solution (m^'^'", z^'^'^) G ^ to the following BSDE: 



u 



N,M,n( 



^N,M,nf^ „N,M,n 



N,M,n/ 



7N,M.n/ 



In a similar way to Lemma 13.11 we deduce that there exists a positive constant Ki 
which is independent of N, M and n such that 



sup \\u 

te[o,r] 



N,M,n 



+ E 



< Ki, a.s. 



(4.10) 



Therefore, taking M > ATi, we have i?M(?x^'*^'"(s)) = 1 and (^i^.A^.«, z^=a^.«) is 
independent of M. Thus, we write (m^'", Z^'"") instead of (u^'*^'", Z^-^^'") below. 
Moreover, there exists a positive constant A'2 independent of n such that 



||$(t,M^'"(t),0i) -$(t,M^'"(t), 02)11 <i^2||0i-02||, rfP ® rft-almost, 
holds for all (f)i,(f)2 ^ H and N,n & Z"*". 



(4.11) 
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For J e Z+, set (f/^,X^) 
similar to (14. 8p . we get 



N,n+j _ ^N,n^ ^A^.n+j _ ^iv,n^_ Applying Ito formula 



2A2 



T 



T 



< 2 / e^^^($^'"+^(s,u^'"(s),Z^'"(s))-<l>^'"(s,M^'"(s),Z^'"(s)), f/^(s))rfs 



-2 j e'^''{U''{s),X''{s))dWs 
(by (moD) 

On the other hand, 



(4.12) 



E 



sup 



< CE 



< eE 



e'''^U''{s),X''{s))dWs\ 



T X 1/2- 

e2^^^||X^(s)|n|f/^(s)f rfs 



(by BDG inequality) (4.13) 



sup (e^^^||f/^(.)f ) 



+ aE 



||X^(s)fe^2^ds 



with the positive constant e to be determined later. Then choosing e to be small 
enough, we deduce from fl4.12p that 

||<I>^'"+^'(S, ^"^'"(5), Z^'"(s)) - $^'"(5, ^"^'"(5), Z^'"(S))|| (is 



As 



(4.14) 



||<I>^'"+^'(S, M^'"(s), Z^'"(s)) - $^'"(S, M^'"(S), ^"^'"(S)) II 

< 2ir2||^'^'"(s)P{,|ZA',n(,)||>„}+2K2||^'^'"(s)P{^,^(.)>„}+2/lKl^ 

in view of (ITTOD and /i/^, e ^^(il x [0, T]), we conclude that (m^'", Z^-") is a Cauchy 
sequence in ^ . Denote the limit by (m^, Z^) G ^ . It is easily checked that 
(m^, Z^) is a solution of the projected problem fl3.2p . 
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Step 3. Estimate (14 .7^ follows from Lemma 13.11 which completes the proof. 

□ 



5 Proof of Theorem 2.1 



Proof of Theorem 12. 1[ Our proof consists of the following four steps. 

Step 1. By Theorem 14. 3[ we have solved the projected problem (13.21) in By 
Lemma [3.11 we have 



ess sup \\u 

{u},s)enx[o,T] 



N, 



\Z'^{s)\\Us 



(5.1) 



< 



C {\\g\\L°°{n,L^{[o,T])) + UWl^inyy) 



where C is a constant depending only on u, \, X, (3, g, pi and T . Since we get 
\\B{v)r < CgW f|| ll'y lly from Lemma [2^21 under Assumptions (Al)-(A3), we 
conclude 



E 

Hence, 



+ ||(a.V)Z^(s)f + 11/(3, n^(.),Z^(s))f) ds 



< C. 



|P^$(-,«^,Z^)|U2(^,[0,T],H) < l|$(-,«^^^) 



\L2{nx[0,T],H)\ 



< C. 



(5.2) 



All the constants Cs above are independent of A^. 

Step 2. Now we consider the weak convergence. Clearly, 



C ^ strongly in V, a.s., and ||^^^ ||y < as A^ cx), 

which implies that — > ^ in Lp{Q,V) for any p G (l,+oo). Then the following 
weak and weak star convergence results in respective spaces hold: there exists a 
subsequence {A'fc}^^ of {A^}, such that, as — )■ cxd. 



U 

<|.(■,«^^z^^)(•) 

P^<|.(-,«^SZ^'=)(-) 



^(■)in^^(0,T;D(A))), 
u{-)iYi LJ(fi;C([0,T];F)), 
Z(-) in L^(fi;L2(0,T; V)), 
r(-)in L%{n-L\Q,T-H)), 
^(■)in L%{n-L\{),T- H)), 
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(5.3) 



where m, Z, T and \1/ are some functions in the respective spaces. 

By the Burkholder-Davis-Gundy (BDG in short) inequahty, we can get 



E 



T 



Z^''{s)dWs 



2 






dt 


< TE 


sup 


V 




te[o,T] 



T 2 

I W.' I i 1 1/1/ 

V 



Z^^{s)dWs 









2 " 








2 " 


< 2TE 


sup 


[ Z^''{s)dWs 




+ 2TE 




[ Z^'={s)dWs 






iG[0,T] 


Jo 


V 






'o 


V 








2 " 








< ATE 


sup 


[ Z^''{s)dWs 




< ALiTE 




ds 




te[o,T] 


Jo 


V 




Jo 





(5.4) 



where Li > is a uniform constant from the BDG inequahty. Hence, as a bounded 

'2 /n. t2i 



hnear operator on the space L^(f2;L (0,T; V)), the mapping 

fT 



T : Z^H 



Z^Hs) dW, 



maps the weakly convergent sequence {Z^ ''{■)} to a weakly convergent sequence 
^j^ Z^^{s)dW,^ in L^(fi;L2(0,T;\/)) such that 



T 



Z{s)dWsm L%{VL]L'^{'d,T]V)), as oo. (5.5) 



Similarly it can be shown that, as — )■ oo. 



^(s)cis in L^((];L2(0,T;if)). (5.6) 



Define 



M(t) =i + ~ ^ ^(^) 



(5.7) 



It is easily checked that u = P ® dt-almost. In view of [13i, Theorem 4.2.5], we 
conclude that u e L°°((], C([0, T],V)) and also. 



ess sup ||M(s)||y + i? 
(w,s)e!rJx[o,T] 



T pT 

2 ^„ , / II v^^Mia 



^0 



+ II^IIl- (n;y) f ' 



< C I ||fi'||L°°(n,Li([0,T])) 

where C is a constant depending on i/. A, A, /3, f), pi and T. 



(5.^ 
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Step 3. For a notational convenience, we now use the index N instead of Nk for 
all the relevant subsequences. 

As U^^i^|(0, T; PnH) is dense in ^Jr(0, T; H) and it can be checked that ^ = T 
on U^^]^=SfJr(0, T; PnH), by a density argument we have = F. Thus, to show that 
{u, Z) is a strong solution of the 2D BSNSE problem (12.41) . we only need prove 



^(■) = Z), a.s.. 



(5.9) 



For any v e L^{n, C([0, T],V)), define 



r{t) =r{u,t) := / (K + -||t;(a;, s) ||^ + Kp2(i;(a;, s))) ds, iu,t) e Q x [0,T], 



where the constant K comes from (I3.14p in Lemma 13.21 Applying Ito formula to 
compute e^(*^||M^(t)|p, we have 



\T)r] 



E 



e^(^) ({2PnHs, u'^is), Z^(s)), «^(.)) - ||Z^(.) 



4, 



-{K + jMs)ry + Kp\vis)))\\u''{s)r) ds 



E 



e^(^) 2($(., u'^is), Z^is)) - $(., t;(.), Z(s)), «^(.) - t;(.)) 



^ + jHs)fy + irp2(^(s)) ) ||m^(s) - vis)e 1 t/s 



||Z^^(s)-Z(s) 



+ E 



^^(^) 2(<|.(., «^(.), Z^(s)) - $(., t;(.), Z(.)), t;(.)) 



+ 2($(.,t;(.),Z(.)),n^^(.))-2(Z^^(.),Z(.)) + ||Z(.)r 
-(K+l\\v{s)ry + Kp\v{s))\ {2{u'^{s)Ms))-Hs)r) ) ds 



< E 



e^W 2($(5, w^(3), Z^(s)) - $(., t;(s), Z(s)), t;(.)) 



+ 2($(.,t;(.),Z(.)),«^^(.))-2(Z^^(.),Z(.)) + ||Z(.)r 

i^+l||<s)f^ + Kp2(^;(.))] {2{u''{s)Ms))-\\v{s)r) )ds 



Letting — )■ oo, by Lemma 13.21 and the lower semi continuity, we have for any 
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nonnegative G L°°(0,T), 



E 



m (e^(*)|k(t)ir-e^(^)l|«(T)f) dt 



< lim inf E 



< E 



(5.10) 



+ 2ms, v{s), Z{s)), u{s)) - 2{Z{s), Z{s)) + \\Z{s) 



^ + + (2M5),^^(5))-||t;(5)f) )ds)dt 



while Ito's formula yields 



E [e'^(*^||M(t)|| 



2 _ pKT) I 



E 



(5.11) 



ds 



By substituting (15.111) into (I5.10p . we get 



E 



^{t) e^(^) (^2(vl/ - $(s, v{s), Z{s)),u{s) - v{s)) 

- (^K + j\\v{s)\\l + Kp'{v{s)))\\u{s)-v{s)\\Ads]dt 



(5.12) 



< 0. 



Take v = u- -f(pw for 7 > 0, w; G 1/ and G L°°(fi x [0, T], ?, R). Then we divide 
by 7 and let 7 — )■ to derive that 



E 



< 0. 



3"("V(s)(2(^ - ^{s,u{s),Z{s)),w)ys^dt 

By the arbitrariness of ip,(f) and w, we have 

r = ^ = $(-, M, Z), a.e. on ^] X [0, T] 
In view of (15. 7p and keeping in mind the fact u = u dt x F-a.e., we have 

Z{s) dWs. 



(5.13) 



u{t) = ^+l ^{s,u{s),Z{s))ds 
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(5.14) 



Hence, by Remark [2] we conclude that {u, Z) is a strong solution to the 2D BSNSE 
problem (12 .4^ . 

Step 4. We shall prove the uniqueness. Suppose that there are two strong 
solutions {u{-), and (w(-), Y{-)) to the problem (12.41) corresponding to the same 
terminal data S,- Then 



u{t) - v{t) 



vAu{s) + vAv{s))ds+ / {B{u{s)) - B{v{s))) ds 



+ {JZ{s)-JY{s))ds+ {f{sMs),Z{s))- f{s,v{s),Y{s)))ds 



{Z{s) ~Y{s))dWs, t e [0,T], a.s.. 



Define 



R{t) = R{u,t)= [ {K + hv{u,s)\\1, + Kp'^{v{u,s)))ds, {u,t) e Q x [0,T]. 
Jo ^ 

Then in view of Lemma 13.21 and by Ito's formula (for instance, see fiSl Theorem 
4.2.5]), we have 

E^Je^W|mt)-t;(t)f] 



--E 



T 



e^W(2($(s,M(s),Z(s))-$(s,t;(s),y(s)),M(s)-t;(s)) 
4 



\\Z{s) - Y{s)r -U+ jMsWv + Kp\v{s)) j ||^(.) - v{s)r) 



ds 



(5.15) 



Thus, 



2A2 



■\\Z{s)-Y{s)r-\\\u{s)-v{s)\\l)ds 



E 



3^W||M(t)-t;(t)||2+ /e^(^) 



2A2 



■\\Z{s)-Y{s)f-,\\\u{s)-v{s)fy\ds 



<0, 



which implies 

for any t G [0, T], M(t) — v{f) = in if, a.s. and E 



\\Z{s)-Y{s)fds 



By the continuity of u and f , we have ||(ti — — y)||_^ = 0, from which we conclude 
that (u, Z) is only a modification of (w, Y) in (L^(l^; C([0, T]; y))n^|(0, T; /^(A))) x 
^1(0, T; V"). We complete the proof. □ 
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